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Two- component CH system: 
Inverse Scattering, Peakons and Geometry 

D. D. Holm 1 and R. I. Ivanov 1 2 



Abstract 

An inverse scattering transform method corresponding to a Riemann-Hilbert problem is 
formulated for CH2, the two-component generalization of the Camassa-Holm (CH) equa- 
tion. As an illustration of the method, the multi - soliton solutions corresponding to the 
reflectionless potentials are constructed in terms of the scattering data for CH2. 
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1 Introduction 

Purpose of the paper: 

In this paper we investigate various aspects of the two-component CH system (CH2), including 
its soliton solutions in the inverse scattering framework. The main difference from the standard 
inverse scattering transform method is that the spectral problem for CH2 is a Schrodinger equation 
with an 'energy dependent' potential that is quadratic in the spectral parameter. 



1.1 CH equation 

This section introduces the CH equation and its two-component extension, CH2. Later sections 
will discuss the isospectral problem for the CH2 system, leading eventually to its multi - soliton 
solutions. 

The CH equation [3, 4] 

u t - u xxt + 2uu x + 3uu x - 2u x u xx - uu xxx = 0, (1.1) 

has gained popularity as an integrable model describing the unidirectional propagation of shallow 
water waves over a flat bottom [3, 4, 12, 13, 31, 32] as well as that of axially symmetric waves in 
a hyperelastic rod [11]. In the shallow water wave interpretation of CH, the real parameter u is 
the asymptotic value of the horizontal fluid velocity u at spatial infinity, as \x\ — > oo. For u = 
the CH equation possesses singular solution in the form of peaked solitons (peakons) [3, 4, 18]. 
Summaries of the developments of many results about the CH equation appear in, e.g., [24, 18, 19] 
and references therein. In the present context its most important properties are its bi-Hamiltonian 
structure and its Lax pair. 

The CH equation in (1.1) may be expressed in bi-Hamiltonian form as 

m t — —[d — a) — — - — - = —(dm + mo + 2cuo) — — -, (1.2) 
dm dm 

where the momentum m associated to the fluid velocity u is given by 

m = u — u xx (1.3) 



and the two Hamiltonians are 



Hi[m] — - I mudx, (1.4) 



2 



H- 



[m] = ^ J (w 3 + uu 2 x + 2uu 2 )dx. (1.5) 
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The integration is taken over the real line for functions rapidly decaying as \x\ — > oo, and taken 
over one period in the periodic case. (In the periodic case, ui is related to the mean depth.) 

The CH equation admits an infinite sequence of conservation laws (multi-Hamiltonian structure) 
H n [m], n = 0, ±1, ±2, . . ., obtainable from the recursion relation 

— (a — a) — — = — (dm + mo + 2uo) . (1.6) 

dm dm 

The recursion relation for CH leads to its Lax pair. Namely, the CH equation follows as the 
compatibility condition for the Lax pair [3, 4] 

*xx = Q + A(m + w))tf, (1.7) 

= (^-«)** + y* + 7*» ( L8 ) 

in which 7 is an arbitrary constant. 



1.2 From CH to CH2 

An integrable two-component generalization of the CH equation can be easily obtained by extend- 
ing the Lax pair for CH in (1.7), (1.8) to a Lax pair whose eigenvalue problem is quadratic in the 
spectral parameter [5] 



*„ = I -aVW + M(x) + - 4 \ *, (1.9) 

*■ = (k -»)**+ t*> 

The compatibility of the two equations (1.9), (1.10) produces a two-component extension of the 
CH equation, abbreviated as CH2, 



qt + uq x + 2qu x + pp x = 0, (1.11) 
Pt + M, = 0, (1.12) 

where q = u — u xx + 00 with 00 being a constant. In our further considerations, we shall assume the 
limit relation lim| a; |_ 5 . 00 (p(x) — po) = 0, where p > is a constant, while both u(x) and p{x) — po 
are Schwartz class functions. Taking p = p = reduces the CH2 system to the CH equation. 

The CH2 energy Hamiltonian is given by 



H =\J [u 2 + u 2 x + {p- p G f]dx, 



(1.13) 
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and is positive-definite. The CH2 system (1.11-1.12) is bi-Hamiltonian. This means it has two 
compatible Poisson brackets. Its first Poisson bracket between two functionals F and G of the 
variables m and p is in semidirect- product Lie-Poisson form [22, 17] 



{F,G}, 



-I 



5F 5G 5F 5G 8F 5G 

- — (md + dm)- h -=—po- — V -r-op- — 

dm dm dm op op dm 



dx. 



;i.i4) 



This Poisson bracket generates the CH2 system from the Hamiltonian H\ 
m = u — u xx . Its second Poisson bracket has constant coefficients, 



{F,G} 2 = 



dF 
dm 



(d-d 3 



dG 
dm 



SF Q 8G 
dp dp 



dx, 



| J (urn + p 2 )dx with 



(1.15) 



and corresponds to the Hamiltonian H 2 — \ J (up 2 + u 3 + uu 2 )dx. There are two Casimirs for the 
second bracket: J pdx and J m dx. Since H and Hi differ only by a Casimir of the second bracket, 
they both generate the same flow (uniform translation: m t + m x = and p t + p x = 0) under the 
Poisson bracket (1.15). 

The CH2 system represents a two-component generalization of the CH equation. It was initially 
introduced in [37] as a tri-Hamiltonian system, and was studied further by others, see, e.g., [35, 5, 
15, 28, 10, 23, 16, 39]. The CH2 model has various applications. For example: 



• In the context of shallow water waves propagating over a flat bottom, u can be interpreted 
as the horizontal fluid velocity and p is the water depth in the first approximation [10, 29]. 

• In Vlasov plasma models, CH2 describes the closure of the kinetic moments of the single- 
particle probability distribution for geodesic motion on the symplectomorphisms [25]. 

• In the large-deformation diffeomorphic approach to image matching, the CH2 equation is 
summoned in a type of matching procedure called metamorphosis [26]. 



For discussions of the geometric aspects of the CH2 system we refer to [24, 34, 26]. Its analytical 
properties such as well-posedness and wave breaking were studied in [14, 10, 27, 40, 16, 6] and 
elsewhere. In general, one can show that small initial data of the CH2 system develop into global 
solutions, while for some initial data wave breaking occurs. Only the plus sign (+) in front of 
the pp x term (1.11) corresponds to a positively defined Hamiltonian and straightforward physical 
applications. It would be interesting to know whether the model with the choice of the minus sign 
in (1.11) has a physical interpretation, since this case is also integrable [10]. 



Solutions of CH2 for dam-break initial conditions. Figure 1.1 plots the evolution of CH2 
solutions for (u, p) governed by equations (1.11-1.12) with the + sign choice in the periodic domain 
[—L,L] with dam-break initial conditions given by 

u(x,0) = 0, p (x, 0) = 1 + tanh(x + a) — tanh(x — a) , (1-16) 

where a<L. 

The dam-break problem involves a body of water of uniform depth, initially retained behind a 
barrier, in this case at x = ±a. When the barrier is suddenly removed at t — 0, the water flows 
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Figure 1.1: Dam-break results for the CH2 system in equations (1.11-1.12) show evolution of the 
density p (left panel) and velocity u (right panel), arising from initial conditions (1.16) in a periodic 
domain. The color bars show positive density on the left and both positive and negative velocity 
on the right. The soliton solutions are seen to emerge after a finite time, and the evolution of 
both variables generates more and more solitons propagating in both directions as time progresses. 
Figures are courtesy of L. O Naraigh. 

downward and outward under gravity. The problem is to find the subsequent flow and determine 
the shape of the free surface. This question is addressed in the context of shallow-water theory, 
e.g., by Acheson [1], and thus serves as a typical hydrodynamic problem of relevance for CH2 
solutions with the + sign choice in (1.11). 

Plan of the paper 

Section 2 discusses the isospectral problem for the CH2 system. Section 3 treats asymptotics of 
the Jost solutions for the CH2 scattering problem. Section 4 explains how analytic solutions for 
CH2 are obtained by formulating the Inverse Scattering Transform (1ST) for CH2 as a Riemann- 
Hilbert problem (RHP). Perhaps not unexpectedly from the viewpoint of CH2 as a fluids system, 
its solutions possess the parameterised form (4.15-4.17) corresponding to fluid continuum flow. 
Section 5 treats multi - soliton solutions of CH2 arising as reflectionless potentials. That is, the 
reflection coefficient in the inverse scattering transform is taken to vanish in the solution of the 
RHP for CH2. Section 6 provides a slightly modified CH2 equation that admits peakon solutions, 
but may not be integrable. Section 7 closes the paper by giving a brief summary of its main points 
and indicating some directions for future research. 

2 The scattering problem for CH2 

Outlook for the CH2 scattering problem. This section begins our discussion of the isospec- 
tral problem for the two-component CH equation with a single velocity, denoted CH2 for simplicity. 
The next three short sections will be devoted to further discussions of the CH2 scattering problem. 
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In §3, we will treat the asymptotics of the Jost solutions for the CH2 scattering problem. In §4, 
we will explain how to formulate the Inverse Scattering Transform for CH2 as a Riemann-Hilbert 
problem. Finally, in §5 we will derive multi - soliton solutions of CH2 that arise as reflectionless 
potentials. 



2.1 CH2 spectral problem 

The spectral problem for CH2 (1.9) is a type of Schrodinger equation with an 'energy depen- 
dent' potential. In particular, it is quadratic in the spectral parameter and the potential functions 
multiply the spectral parameter. (This is the so-called weighted problem.) It shares some fea- 
tures in common with Sturm-Liouville spectral problems, see for example [30, 33, 38]. An 'energy 
dependent' spectral problem also appears in the inverse scattering transform of an integrable gen- 
eralization of the Bousinesq equation (Kaup-Bousinesq equation) [33]. 

Asymptotically, as \x\ — > oo, the spectral problem (1.9) for CH2 reduces to 

*xz = [-P 2 \ 2 + oj\+^, (2.1) 

or, simply 

*xx = -k 2 *, (2.2) 
where we introduce a spectral parameter k via the equation 



-pl\ 2 + u\ + + k 2 = 0, (2.3) 

The solutions of (2.2) oscillate for real k. Consequently, the continuous spectrum is the real line 
in the complex /c-plane. 

The quadratic equation (2.3) has roots, 



us ok I Pn + co 2 

A <^=^ + W 1+ w- i2A) 



where a = ±1, and we assume t^/w = ^\w\e2 Ar ' 9< - wS) where < Arg{w) < 2ir. An expansion of 
(2.4) for large \k\ yields, 

Mk '° )= «, + M + ^rk +0 (>?y (2 - 5) 

This expansion uniquely determines A from k and a. Equation (2.4) possesses a reflection property 
that we will assume explicitly from here on, that 

X(-k,-a) = X(k,a), (2.6) 
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and also, for real k, 

X(k,a) = X(k,a), (2.7) 

where A is the complex conjugate. 

As usual, for real k ^ a basis in the space of solutions of (1.9) can be introduced, fixed by its 
asymptotic behavior when x — > oo [36, 7, 8]: 

il>i(x, k) = e- ikx + o{l), x^oo; (2.8) 
^ 2 (x,k) = e ikx + o{l), x-^oc. (2.9) 

A complementary basis can also be introduced, fixed by its asymptotic behavior when x — > — oo: 

ifi(x,k) = e- ikx + o(l), x-+-oo; (2.10) 
<p 2 (x,k) = e ikx + o(l), x^-oo. (2.11) 

2.2 The scattering matrix, the Jost solutions and the reflection coef- 
ficient 

Since A depends not only on k but also on a, it follows that the entire spectral problem, as well 
as the eigenfunctions, are labelled by a. For all real k ^ 0, if ^/(x, k, a) is a solution of (1.9), then 
fy(x, —k, —a) is also a solution, since they share the same A, according to (2.6). Thus, 

(pi(x, k, a) — <pi(x, —k, —a), tpi(x, k, a) = vp 2 (x, -k, -a). (2.12) 
Due to the reality of q, p in (1.9) and the property (2.7) for A, for real k we have 

(fi(x,k,a) = Lp 2 (x,k,a), ipi(x, k, a) = 4> 2 (x, k, a). (2.13) 

For real k the vectors of each of the bases may be represented as a linear combination of the vectors 
of the other basis: 

<Pi(x, k,a) = ^ T u( k i a )M x ^ k , °")> (2-14) 

1=1,2 

where the matrix T(k,a) defined above is called the scattering matrix. 

For real k ^ 0, instead of <fi(x, k, a) and (p 2 (x,k,a), due to (2.13), for simplicity we can write 
correspondingly ip(x, k, a), {p(x,k,a). Similarly, we can replace ipi t2 (x,k,a) by ip(x,k,a) and its 
complex conjugate ^(x, k, a). Thus T(k, a) has the form (with real k) 
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and clearly 

ip(x, k, a) = a(k, a)ip(x, k, a) + b(k, a)ip(x, k, a). (2-16) 

Remark. The solutions ip(x, k, a) and ip(x, k, a) are called the Jost solutions. 

The Wronskian W(fi, / 2 ) = fid x f 2 — /^/i of any pair of solutions of (1.9) does not depend on 
x. Therefore, perhaps not unexpectedly, 

W((p(x, k, a), <p(x, k, a)) = W(i/j(x, k, a),^(x, k, a)) = 2ik. (2.17) 
From (2.16) and (2.17) it follows that 

\a(k,a)\ 2 - \b(k,cr)\ 2 = I. (2.18) 

Hence, det(T(k,a)) = 1. That is, the determinant of the scattering matrix is unity. 

In analogy with the spectral problem for the KdV equation, (which is the Schrodinger equation 
from quantum mechanics) [36], one can introduce a reflection coefficient, 

TZ(k, a) = b(k, a)/a(k, a). (2.19) 

The matrix T(k,a) in (2.15) is determined from the knowledge of 1Z(k,a) for real k > only. 
Indeed, from (2.12) we have ip(x, k,a) = tp(x, —k, —a), etc. for real k. Hence, the scattering data 
satisfy 

a{k,a) = a(-k,-a), b(k,a) = b(-k, -a). (2.20) 
Also, it is sufficient to know TZ(k, a) only on the half line k > 0, since 1Z(—k, —a) = TZ(k, a). 

3 Asymptotics of the Jost solutions for CH2 as \k\ — > oo 

Outlook. This section continues the analysis of the CH2 scattering problem, by discussing the 
asymptotic behavior of the Jost solutions. 

3.1 Analyticity properties 

The analyticity properties of the Jost solutions and of a(k, a) play an important role in our con- 
siderations. We will need also the asymptotic behavior of the Jost solutions for \k\ — > oo which 
have the form (cf. [8, 9]) 



X (x) + O 



(3.1) 
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ip(x, k,a) — e 



X (x) + O 



k 



(3.2) 



where Xn(x) 



1/2 



. (If initially p(x, 0) > 0, one can easily prove that p(x,t) always remains 

positive.) As in [8] one can show that ip(x,k,cr) is analytic in the lower complex half k -plane, 
while <p(x, k,a) is analytic in the upper complex half k -plane. 



The expression for a(k, a) can be extended into the upper half plane by 

a(k, a) = ^ W (f( x ^ k , <t),iJ>(x, ~k, -a)). 
An immediate consequence of (3.1), (3.2) and (3.3) is: 



(3.3) 



where the quantities 



lim a{k,ay ka - i(lfi = 1, k G C + , 

k— >oo 



a = 



p(x) 

-oo \ PO 



1 dx, 



q(x) up(x) 



dx. 



(3.4) 

(3.5) 
(3.6) 



are two integrals of the system, 



3.2 The discrete spectrum. 

The discrete spectrum can be found as follows. Suppose that ko(a) G C + is a zero of a(k, a). Then 
ip(x, k , a) and ip(x, —ko, —cr) are linearly dependent (2.16): 

ip(x, k , a) = b (a)tp(x, -k , -a). (3.7) 

From here we see that ip(x,ko,&) decays exponentially for both x — > — oo (which follows from 
the definition of ip(x, k , cr)) and x — >■ +oo (since ip(x, — k , — a) = e lkx for x — > oo). Therefore 
ip(x, k , a) is a well defined eigenfunction of the discrete spectrum with an eigenvalue ko- 

Now, multiplying (1.9) by (p{x, ko,a) and performing some manipulations while keeping in mind 
that the eigenfunction decays exponentially for both x — > ±oo, we obtain 



/OO POO 
q 2 (x)\if\ 2 dx + A(fco) / qi(x)\if\ 2 dx + 
-oo J — oo 

j (hipi 2 + \i Px \ 2 ')dx = o 

J —oo ^ 



(3.8) 
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This identity can be regarded as an equation for X{ko) where ko is a parameter. From the quadratic 
formula, the two roots X(ko, cr) and X(ko, —cr) satisfy 

X(k , a)X(k , -a) = ^ 2( r, ' ■ (3.9) 

Loop \ x )m Ax 

On the other hand, from (2.3) we have 



X(k ,a)X(k ,-a) = -\fk 2 + ]). (3.10) 



From (3.9) and (3.10) it follows that k^ + \ is real and positive. Since ko is in the upper half complex 
plane, it should be exactly on the imaginary axis, ko = i^o, where k is real, and < k < 1/2. 
With this restriction on ko, notice that X(ko,a) is real. 

Let us show that a(k, a) can have only simple zeroes in the upper half complex plane. The dot 
will be used to denote the derivatives with respect to k at the point k . From (2.3) we have 
A = 2/c/(2pqA — to). Differentiating the (1.9) (written for the eigenfunction ip) with respect to k 
and multiplying by (p we obtain 



/oo 
(g-2Ap 2 )M 2 dx. (3.11) 
■oo 

Next, using the asymptotics ip(x,ko,cr) — > a(k ,a)e K ° x , ip(x,ko,cr) — > b {a)e^ K ° x for x — > oo; 
ip(x, k , a) — > 0, ip(x, k , cr) — > for x — > — oo, (3.11) can be transformed into 

/oo 
(q - 2A(A: )p 2 )|^| 2 drc = (u - 2p 2 X(ko))ibod(k ). (3.12) 
-oo 

As a corollary we notice that the quantity Ro(cr,t) = -r^y is real. 

If f^dq — 2X(ko)p 2 )\(f\ 2 dx ^ 0, then a(k ) ^ and the zero k is simple. Therefore, a multiple 
zero is possible, only if 

/oo 
(g-2A(A; )p 2 )M 2 d:r = 0. (3.13) 
-oo 

Suppose that (3.13) is satisfied. From (1.9) (written for the eigenfunction <p) multiplied by <p we 
obtain 

/oo poo 
p 2 (x)\<p\ 2 dx + X{k ) / q(xM 2 dx 
■oo J — oo 

+ [ QM 2 + |y^| 2 W = 0. (3.14) 

J —oo ^ 
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From (3.13) and (3.14) we obtain 



Kh) 2 = 



J-oo ^4 



4 iV?| 2 + \¥x\ 2 ^dx 



I- 0O P 2 (x)\f\ 2 dx 



(3.15) 



but from (3.14) itself we find that the product of the two roots is 



rJiM 2 + K| 2 )d:r 
A A*, <r A fc„, -a) = ^ 2 , | 2 / . (3.16) 

From (3.15), (3.16) and (2.4) we find that X(k ,a) = X(ko,—a) = Thus in this case the 

multiplier u — 2pqA(A; ) on the right hand side of (3.12) is also zero. Therefore, we can use 
l'Hospital's rule in the evaluation of a(ko). From (3.12) we have 

AT ( 

d(fco) = lim dxJ -°° 



since 



9 JZia ~ 2\p 2 M 2 dx 
aT 5j fb £(uj-2\p 2 ) 
Um -2 JZ P 2 M 2 dx + 2g jrjg - 2\p 2 M Ijgjdx 

-2ib p 2 

iboPo 



dk 2p 2 \ - u 
lim — = lim = 0. 

dX A— 2k 

2 P0 2 "0 



(3.17) 



Then (3.17) shows that a(k ) ^ 0, i.e. k is a simple zero of a(k). 



3.3 Summary of asymptotic behavior of Jost functions for CH2 

To summarize: the discrete spectrum in the upper half plane consists of finitely many points 
k n = iK n , n = 1, . . . , N, which are the simple zeroes of a(k,a). Furthermore, each K n is real and 
< n n < 1/2. 

Eigenfunctions. Two eigenfunctions <p( n \x, a) belong to each eigenvalue m n , because there are 
two eigenvalues A n (<r) = X(iK n ,a) that correspond to a given K n . We can take this eigenfunction 
to be 

ip( n \x, a) = ip(x, in n , a) . (3.18) 
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The asymptotic behavior of (p( n \ according to (2.10), (2.9), (3.7) is 

ip^(x,a) = e KnX + o(e KnX ) for x ->■ -oo , (3.19) 
^ n \x,a) = b n (a)e~ K " x + o(e- KnX ) for x ->■ oo. (3.20) 

Scattering data. The set 

5 = {7£(A;, cr) (A; > 0), K n , b n (a), n = l,...N, a = ±1}, (3.21) 

is called the scattering data. 

The time evolution of the scattering data can be easily obtained as follows. 
The second equation of the Lax pair is 

*t = (^-«)** + y* + 7*> (3-22) 

where we introduced an arbitrary constant 7 (which does not affect the compatibility). 
From (2.16) with x — > 00 one has 

ip(x, k, a) = a{k, a)e~ ikx + b{k, a)e ikx + o(l). (3.23) 

The substitution of A;, cr) into (3.22) with x — > 00 gives 

¥>t = 7^<Ac + 7V? • (3-24) 
From (3.23), (3.24) with the choice 7 = ik/2\ for the eigenfunction ip(x, k) we obtain 



a t (k,a,t) 
b t (k,t,a) 

Thus, we find 

a(k, t, cr) = a(k, 0, cr), 



= 0, (3.25) 

= ^b(k,t,a). (3.26) 
A 

6(Jfe, cr) = b(k,0,a)e^ t , (3.27) 

(3.28) 



TZ(k, t, cr) = 7£(fc, 0, jJeMW 4 . 



In other words, a(/c, cr) does not depend on t and can serve as a generating function of the conser- 
vation laws. 
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Time evolution of the data on the discrete spectrum. The time evolution of the data 
on the discrete spectrum is found as follows. Let us introduce the notation a n (a) = a(iK n ,a). 
We notice that iK n are zeroes of a(k,a), which does not depend on t, and hence (n n )t = and 
Ki{&)t — 0. From (3.22), (2.4) with 7 = ik/2\; k = ik u and (3.20) one can obtain 

b n (a) t = ^\b n {a). (3.29) 

It is convenient to use the variable R n (cr) = j^j^j, which is real, according to (3.12) and evolves 
with t as 

R n (t, a) = R n (0, a) exp (^rt) . (3.30) 

When k is in the upper half plane one can derive the following dispersion relation for a(k, a) e.g. 
following the pattern for the CH case from [9]: 

m*. *)=-*•* ^j^+hC^HT*' (3 - 31) 

i.e. a(k, a) is determined by \a(k, a)\ given on the real line k e M. 



Outlook. In the next section, we will develop the Inverse Scattering Transform for CH2. The 
special case of reflectionless potentials (lZ(k,a) = for all k) corresponds to an important class of 
solutions, namely the multi-soliton solutions. These will be separately studied in Section 5, where 
a formula for the iV-soliton solution will be obtained. 



4 Analytic solutions and the Riemann-Hilbert Problem for 
CH2 

This section explains how analytic solutions for CH2 are obtained by formulating the Inverse 
Scattering Transform for CH2 as a Riemann-Hilbert problem (RHP). 

4.1 Preliminaries 

We begin by introducing the following new variables, cf. the integrals of motion in equation (3.5) 
and (3.6) 

y{x) = x+ I* (^-l)dx', (4.1) 
In terms of these variables, the expansion (3.1) may be written as 
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ip(x,k,a) = e - lky+iaz 



X (x) + O 



(4.3) 



Furthermore, the function x( x ,k,a) = ■ifj(x,k,o-)e lkx is analytic for Im k < 0, due to arguments 
similar to these, given in [7] for the CH case. This follows from the representation 



/oo 2ik(x—x') i 
\-X\p\x') - pi) + \{q{x>) - u))\ x(x', k)dx>. 



2ik 



Notice that y(x) is a bounded function for all x, which follows from the assumption that p{x) — p Q 
is a Schwartz class function. Therefore, the function 



f(x, k, a) = i>(x, k, a)e iky ~ iaz = X (x) + O (^j 



is also analytic for Im k < 0. 
Similarly, 



(p(x, k, a) = <p(x, k, a ) e iky-i*z+ika-i*f3 = Xo ( x j + q 



(4.4) 



(4.5) 



is analytic for Im k > 0. 

Multiplying (2.16) by e iky ~ icTZ /a(k, a) and using (4.4) and (4.5) we obtain 



ri! ' / '"' i/>(x,k,a) +n{k,o-)^{x,-k,-(r)e 2iky - 2icTZ . (4.6) 



e ika - ia Pa(k, a) 



The function eik ^^,'^ k ^ is analytic for Im k > 0, while ^{x, k, a) is analytic for Im k < 0. Thus, 
equation (4.6) represents an additive Riemann-Hilbert Problem (RHP) with a jump on the real line, 
given by 1Z(k, o-)ijj(x, —k, —a)e 2lky ~ 2mz and a normalization condition, limifc^oo ip(x, k, a) = X (x). 



4.2 Solving the Riemann-Hilbert Problem for CH2 

In this section we will follow the standard technique for solving RHP. We integrate the two an- 
alytic functions with respect to § ^z^( - ) over the boundary of their analyticity domains, using 
the normalization condition. In our case the domains (the upper C + and the lower C_ complex 
half-planes) have the real line as a common boundary and there we relate the integrals using 
the jump condition. The RHP approach for the CH equation is presented in [8, 2, 19], for the 
Kaup-Bousinesq equation (which also has energy-dependent spectal problem) in [38] . 
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Let us take an arbitrary k from the lower half plane (Im k < 0). Then using the Residue Theorem, 
(3.5) and (3.7) we can compute the integral 



I = 



tp(x, k', a) dk' 
2m J c+ e ik ' a - ia Pa(k> ', a) k' - k 



N 

E 



tp(x,iK n ,a) 



^ (iK n - k)e- Kna - iaf} a n (a) 
^ iR n (a)e- 2K -y- 2i ^{x, -iK n , -a) 



n=l 



(4.7) 



where C + is the closed contour in the upper half plane (Fig. 4.1). 



+ 



© 



Figure 4.1: The contours r± 



On the other hand, because of (4.6) the same integral can be computed directly as 



l - J (±(x, k', a) + U(k', a)±(x, -k', -o)e™'^° z ) ^ 
<f(x, k', a) dk' 



+ 



2ni 



2?ri J T e ik ' a ~ i ' J f i a{k' \a) k' - k' 



(4.8) 
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where T + is the infinite semicircle in the upper half plane (Fig. 4.1). Using the expansion (4.5) 
and limit (3.4), it is straightforward to compute that the integral over T + is simply (l/2)X (x). 

Similarly, 

- t (x,k,a) = ^J c t^k\*)^ = ^J^±(x,k',a)-? k 



k'-k 



1 I" rlk' 

where C_ is the closed contour in the lower half plane, T_ is the infinite semicircle in the lower 
half plane (Fig. 4.1). Due to (4.3), (4.4) the integral over T_ is -(l/2)X (x). 

Now, from (4.7) - (4.9) it follows that for Im k < 0, 

/°° rlk' 
TZ(k', <j)±(x, —k', -ay 2 **-*"*-— 
■oo k 



k 

(x, -iK n , -a) 

The expression (4.10), taken at k — —iK p , p — 1, . . . , N gives 



" tR n (a)e' 2 ^y- 2 ^l(x, -iKn, -a) yX 

n=l 



i/)(x,—iK p ,cr) = X (x) 

/°° rlk' 
K(k',o-)l(x,-k',-a)e 2ik ' y - 2iaz 
-oo 



k' + 

^ R n (a)e- 2K "y~ 2 ^(x, -m n , -a) ^ ^ 

n=l 1 



From (4.11) it follows immediately that 



^_(x,-iK p ,-a) = X (x) 

/oo 
K(k',-cr)ip(x,-k' 
■oo 



rlk' 

\ c -2ik'y+2iaz aK 



k' + 

^ R n (-a)e - 2K ^ +2i ^iP(x, -iK n , a) 



(4.12) 



Equations (4.10) - (4.12) represent a linear system, from which ip(x, k, a) (for real k) and ip(x, —in n , ±0") 
can be expressed through y, z, which are as yet unknown functions of x (since X (x) can be ob- 
tained from y(x)). 
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Finally, we need to find the dependence of y on x. From the quadratic roots in (2.4) we notice 
that 

A(-i/2,<t)= 1 ' . 

z Po 

Hence, if we take a = <7i = sign (a;), then we have A(— i/2,a\) = 0. Now ip(x,k,ai) does not 
depend on q and p for A = and since ip(x, k, <j\) is defined by its asymptotic behavior as x — > oo, 
which is e~ x l 2 , i.e., it is real when k — — i/2. Consequently, we have 

- (ip(x, —i/2, (Ti) + compl. conj.) = e~ x ^ 2 . 
2 

Thus, for k = —i/2, a — u\ = sign(wi), equation (4.10) gives 



e (-x+y)/2 _ e ia ± z 



(X {x) + J TZ(k', aJKx, -k! , -o x )e™ 'v-^-^-^ 

R n (o-i)e- 2K ™y- 2i ^ z ip(x, -in n , -<ri) \ 
— -|_~ J ~*~ com Pl ex conjugate (4.13) 

In other words, (4.10) - (4.13) represent a system of singular integral equations for ip(x, k, a) (for 
real k) and ip(x, —iK n {a), a), ip(x, —iK n (a), a). 

A similar relation to (4.10) can be written if k e C + and in particular when k = i/2 we have 
another equation of the type of (4.13). Thus, one can recover 

VbX 

x = X(y,z(y)) and — = X 2 (y, z(y)). (4.14) 
Then eliminating z eventually yields x — X(y). 

Parametric forms of the solutions. Since the time evolution of the scattering data is known 
(3.30), the dependence on t, i.e. x = X(y,t), is also known, as expressed by the scattering data. 
Thus the set S of scattering data (3.21) uniquely determines the solution. From (4.1) one obtains 
the parametric forms of the solutions, 

x = X(y,t), (4.15) 

P(M) = Y^—y (4.16) 

Xy{y,t) 

u{x,t) = X t (y,t). (4.17) 

These correspond precisely to the relations between Eulerian and Lagrangian variables in com- 
pressible ideal fluid dynamics. 



5 Reflect ionless potentials and solitons for CH2 

This section discusses the simplification of the inverse scattering problem for CH2 in the important 
case of the so-called reflectionless potentials, when the scattering data is confined to the case 
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lZ(k, a) — for all real k. This class of potentials corresponds to the multi - soliton solutions of 
the two component CH equation. 

The time evolution of R n due to (3.30) is 



R n (t,a) = R n (0,a)ex P [y^)*) " (5 - X) 
As we already observed, both X n (a) and R n (t, a) are real. Let us define the N x N matrix 

M pq (y,t,a) = 5 pq -^ , (5.2) 

n=l ^ p ~r ^1) 

which is real, and the vectors (with real components) 

N 

A n (y,t,a) = ^[M-^cx)]^, (5.3) 
p=i 

B n ( y ,t,a) = Y^[M-\ y ^)\J^£ — ! - ( 5 - 4 ) 



p=l g=l 

The solution of the system (4.11), (4.12) is 

±{x,-iK n ,a) = X (x)[A n (y,t,a)-e- 2laz B n (y,t,a)], (5.5) 

and from (4.10) also 



±(x,k,a) = X (x)[l-J2 



V 2 J R n (t,a)e- 2K ^ J B n ( 2 /,t,-a) 



n=l 

-2j(72 



, „y^ ^n(t,a)e- 2 ^A t (y,t,-g) i 



n=l 



n 



From (4.13) we have 



Rnio-Je-^iA^y, t, -a,) - B n (y, t, -a,))- 



e (-x+y)/2 = cos{aiz)Xo (x)(l-J2 

(5.7) 



^n + 1/2 



A similar relation for ip(x, k, a) can be written if k G C + and for k — i/2 it gives 



Holm & Ivanov 



Multi- component CH equations 



19 



N 



,(*-v)/2 



= cos((7iz)X (x)ll-2^ _ 1/9 ), 

n=l n / 



(5.J 



and finally from (5.7) and (5.8): 



x = X(y,t)=y + \n I j-^ y (5.9) 



n=l 



The time evolution of R n is known (5.1) and thus X(y,t) is given in terms of the scattering data. 
This produces a parametric representation of the solution in terms of the dependent variables, 
from equations (4.15) - (4.17). 

For the one-soliton solution, equations (5.9) and (5.10) yield the parametric (Lagrangian) repre- 
sentation of the "particle paths" , 

1 - 1 p^ E (a 1 )E(-a 1 ) + ^(-n) 
X(y, t) = y + \n ^ ^ , (5.11) 

where E{a,t) = ^- i R 1 (t,a)e- 2Kiy . 

We notice that in (4.13) one can take the imaginary part of the eigenfunction which is also equal 
to e~ x / 2 , up to multiplicative constant (since both the real and imaginary parts decay to zero as 
x — > oo and when A = they both equal e~ x / 2 ). The imaginary part yields another solution, which 
differs from the presented one by change of the signes of the scattering data Ri(t,a). 

Thus, the Inverse Scattering Transform method yields the parametric representation (5.9) - (5.10) 
of the solution of the CH2 system in terms of the scattering data for its isospectral problem. From 
the viewpoint of fluid dynamics, this is a Lagrangian (particle) representation of the solution, 
which may be written in terms of the Eulerian (spatial) representation by using equations (4.15) 
- (4.17). ' 



6 A modified version of CH2 arising from the 2D EPDiff 
equation 

In this section we discuss a modified version of CH2 that does not lie within its integrable hierar- 
chy but does admit peakon solutions. Namely, we consider solutions of the 2-component EPDiff 
equation that depend only on the first spatial variable x = X\ and do not depend on x 2 : 

qt + uq x + 2qu x + p(l - d 2 x )~ ] l p x = 0, (6.1) 
Pt + (up) x = 0. (6.2) 
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This system of equations has been considered previously [23] from another viewpoint and it is 
known that it has peakon solutions. However, a totally different interpretation of these solutions 
exists. Instead of two different types of variables u and p :— (1 — d 2 .)~ 1 p, one a velocity and one 
an average density, we may imagine having two velocity components U\ = u and u 2 = p. With 
this interpretation, the modified CH2 equations (6.1), (6. 2) in ID are equivalent to the original 
EPDiff( J ff 1 ) equation in 2D coordinates (x u x 2 ) E R 2 , [21] 

m t + u ■ Vm + (Vu) T • m = , (6.3) 

with 2D momentum m = u — Au = (mi,rri2) T and velocity u = (ui,U2) T independent of the 
second coordinate, x 2 . Substitution of m(xi,t) and u(xi,t) into the 2D EPDiff(i/ 1 ) equation (6.3) 
yields the ID two-component system (6.1), (6. 2) with x = X\. 

The presence of peakon solutions of (6.1), (6. 2) is not a surprise, since the singular solutions are 
a characteristic feature for the EPDiff equation, which is not known to be integrable beyond its 
one-dimensional version which coincides with the CH equation. The peakon interactions can be 
studied numerically. An example is presented in Figure 6.1, which shows a numerical simulation 
of the velocity u 2 (x\,t) and the particle path x\{t) for a single Lagrangian fluid parcel obtained 
from dxi/dt = Ui(xi(t),t), obtained as a solution of the EPDiff(if 1 ) equation (6.3) in which the 
solutions (ui,u 2 ) are independent of the second coordinate. 




Figure 6.1: The figure shows the particle path x\{t) for a single Lagrangian fluid parcel obtained 
from dx\jdt = Ui(xi(t),t), upon identifying the modified CH2 solutions (u,p) of the ID dam- 
break experiments studied in [23] with solutions (ui,u 2 ) of the EPDiff^ 1 ) equation (6.3) that 
are independent of the second coordinate, x 2 , in a periodic domain. A fluid parcel initially offset 
to the right of center in x\ accelerates gently to the right as the dam-break produces pulses that 
propagate away on the right side. Later, when the leading pulse returns leftward, the parcel is 
pushed suddenly back with considerably greater leftward acceleration and ends up nearly where it 
started. Its trajectory shows that the wave pulses propagate relative to the Lagrangian parcels, so 
they are not frozen into the fluid motion. However, the influence of the returning pulse shows that 
the Lagrangian motion is still strongly coupled to the waves. Figure courtesy of L. O Naraigh. 
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7 Conclusions/Discussion 

Main results of the paper. The stage was set in Section 2 for our discussions of the inverse 
scattering transform method for the two-component CH2 system. The CH2 Jost solutions were 
obtained in Section §3 and their asymptotic behavior was usedin Section §4 to reformulate the 
scattering problem as a Riemann-Hilbert problem (RHP) . By solving the RHP, multi - soliton 
solutions of CH2 were obtained as reflectionless potentials in Section §5. The soliton solutions 
of CH2 arising from the RHP expressed themselves in a parametric form corresponding to the 
Lagrangian representation of fluid dynamics. A slightly modified version of CH2 was found in 
Section §6 by considering translation invariant EPDiff solutions. The peakon solutions of EPDiff 
(diffeons) were shown graphically to interact with the Lagrangian fluid parcels by briefly sweeping 
them along the peakon trajectory. It remains an open problem, as to whether the modified version 
of CH2 in (6.1) - (6.2) is integrable. 

Explicit outstanding problems. The integrable systems properties CH2 discussed here open 
the door for further generalizations and applications, some of which have been presented in the 
paper and others that will be discussed elsewhere. In particular, CH2 is a member of a large family 
of integrable multi-component PDE based on the Schrodinger equation with an energy-dependent 
potential, as discussed in [20]. The numerical simulation of these PDE, and the formulation and 
analysis of their discrete integrable versions can be expected to attract considerable attention in 
future endeavors. One may expect the continuing interest in wave-breaking analysis for CH and 
CH2 to apply to the integrable PDE in the rest of the CH2 family, as well. 
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